Toward clinical exomes in diagnostics and management of male infertility.

Infertility, affecting ∼10% of men, is predominantly caused by primary spermatogenic failure (SPGF). We screened likely pathogenic and pathogenic (LP/P) variants in 638 candidate genes for male infertility in 521 individuals presenting idiopathic SPGF and 323 normozoospermic men in the ESTAND cohort. Molecular diagnosis was reached for 64 men with SPGF (12%), with findings in 39 genes (6%). The yield did not differ significantly between the subgroups with azoospermia (20/185, 11%), oligozoospermia (18/181, 10%), and primary cryptorchidism with SPGF (26/155, 17%). Notably, 19 of 64 LP/P variants (30%) identified in 28 subjects represented recurrent findings in this study and/or with other male infertility cohorts. NR5A1 was the most frequently affected gene, with seven LP/P variants in six SPGF-affected men and two normozoospermic men. The link to SPGF was validated for recently proposed candidate genes ACTRT1, ASZ1, GLUD2, GREB1L, LEO1, RBM5, ROS1, and TGIF2LY. Heterozygous truncating variants in BNC1, reported in female infertility, emerged as plausible causes of severe oligozoospermia. Data suggested that several infertile men may present congenital conditions with less pronounced or pleiotropic phenotypes affecting the development and function of the reproductive system. Genes regulating the hypothalamic-pituitary-gonadal axis were affected in >30% of subjects with LP/P variants. Six individuals had more than one LP/P variant, including five with two findings from the gene panel. A 4-fold increased prevalence of cancer was observed in men with genetic infertility compared to the general male population (8% vs. 2%; p = 4.4 × 10-3). Expanding genetic testing in andrology will contribute to the multidisciplinary management of SPGF.

The expression profiles of cyp19a1, sf-1, esrs and gths in the brain-pituitary during gonadal sex differentiation in juvenile Japanese eels.

Eels are gonochoristic species whose gonadal differentiation initiates at the yellow eel stage and is influenced by environmental factors. We revealed some sex-related genes were sex dimorphically expressed in gonads during gonadal sex differentiation of Japanese eel (Anguilla japonica); however, the expression of sex-related genes in the brain-pituitary during gonadal sex differentiation in eels is still unclear. This study aimed to investigate the sex-related gene expressions in the brain-pituitary and tried to clarify their roles in the brain and gonads during gonadal sex differentiation. Based on our previous histological study, the control eels developed as males, and estradiol-17ß (E2) was used for feminization. Our results showed that during testicular differentiation, the brain cyp19a1 transcripts and aromatase proteins were increased significantly; moreover, the cyp19a1, sf-1, foxl2s, and esrs (except gperb) transcripts in the midbrain/pituitary also were increased significantly. Forebrain gnrh1 transcripts increased slightly during gonadal differentiation of both sexes, but the gnrhr1b and gnrhr2 transcripts in the midbrain/pituitary were stable during gonadal differentiation. The expression levels of gths and gh in the midbrain/pituitary were significantly increased during testicular differentiation and were much higher in males than in E2-feminized females. These results implied that endogenous estrogens might play essential roles in the brain/pituitary during testicular differentiation, sf-1, foxl2s, and esrs may have roles in cyp19a1 regulation in the midbrain/pituitary of Japanese eels. For the GnRH-GTH axis, gths, especially fshb, may be regulated by esrs and involved in regulating testicular differentiation and development in Japanese eels.

A conserved NR5A1-responsive enhancer regulates SRY in testis-determination.

The Y-linked SRY gene initiates mammalian testis-determination. However, how the expression of SRY is regulated remains elusive. Here, we demonstrate that a conserved steroidogenic factor-1 (SF-1)/NR5A1 binding enhancer is required for appropriate SRY expression to initiate testis-determination in humans. Comparative sequence analysis of SRY 5' regions in mammals identified an evolutionary conserved SF-1/NR5A1-binding motif within a 250 bp region of open chromatin located 5 kilobases upstream of the SRY transcription start site. Genomic analysis of 46,XY individuals with disrupted testis-determination, including a large multigenerational family, identified unique single-base substitutions of highly conserved residues within the SF-1/NR5A1-binding element. In silico modelling and in vitro assays demonstrate the enhancer properties of the NR5A1 motif. Deletion of this hemizygous element by genome-editing, in a novel in vitro cellular model recapitulating human Sertoli cell formation, resulted in a significant reduction in expression of SRY. Therefore, human NR5A1 acts as a regulatory switch between testis and ovary development by upregulating SRY expression, a role that may predate the eutherian radiation. We show that disruption of an enhancer can phenocopy variants in the coding regions of SRY that cause human testis dysgenesis. Since disease causing variants in enhancers are currently rare, the regulation of gene expression in testis-determination offers a paradigm to define enhancer activity in a key developmental process.

NR5A1 and NR5A2 regulate follicle development in chicken (Gallus gallus) by altering proliferation, apoptosis, and steroid hormone synthesis of granulosa cells.

Chicken ovarian follicle development is regulated by complex and dynamic gene expression. Nuclear receptor 5A1 and 5A2 (NR5A1 and NR5A2, respectively) are key genes that regulate steroid hormone production and gonadal development in mammals; however, studies on follicular development in the chicken ovary are scarce. In this study, we investigated the functions of NR5A1 and NR5A2 on follicle development in chickens. The results showed that the expression of NR5A1 and NR5A2 was significantly higher in small yellow follicles and F5. Furthermore, the expression of NR5A1 and NR5A2 was significantly higher in follicular tissues of peak-laying hens (30 wk) than in follicular tissues of late-laying hens (60 wk), with high expression abundance in granulosa cells (GC). The overexpression of NR5A1 and NR5A2 significantly promoted proliferation and inhibited apoptosis of cultured GC; upregulated STAR, CYP11A1, and CYP19A1 expression and estradiol (E2) and progesterone (P4) synthesis in GC from preovulatory follicles (po-GC); and increased STAR, CYP11A1, and CYP19A1 promoter activities. In addition, follicle-stimulating hormone treatment significantly upregulated NR5A1 and NR5A2 expression in po-GC and significantly promoted FSHR, CYP11A1, and HSD3B1 expression in GC from pre-hierarchical follicles and po-GC. The core promoter region of NR5A1 was identified at the -1,095- to -483-bp and -2,054- to -1,536-bp regions from the translation start site (+1), and the core promoter region of NR5A2 was at -998 to -489 bp. Two single nucleotide polymorphisms (SNP) were identified in the core promoter region of the NR5A1 gene, which differed between high- and low-yielding chicken groups. Our study suggested that NR5A1 and NR5A2 promoted chicken follicle development by promoting GC proliferation and E2 and P4 hormone synthesis and inhibiting apoptosis. Moreover, we identified the promoter core region or functional site that regulates NR5A1 and NR5A2 expression.

[Analysis of a child with 46,XY Disorder of sex development due to a novel variant of NR5A1 gene].

OBJECTIVE: To analyze the clinical features and genetic basis of a child with Disorder of sex development (DSD). METHODS: A child who was admitted to the Linyi People's Hospital for primary amenorrhoea on July 29, 2019 was selected as the study subject. Clinical data of the child was collected. Chromosomal karyotyping and quantitative real-time PCR were used to detect Y chromosome microdeletions and other chromosomal aberrations. Next-generation sequencing was carried out for the child and her parents. Candidate variant was verified by Sanger sequencing and bioinformatic analysis. RESULTS: The child, a 13-year-old girl, has featured primary amenorrhoea and onset of secondary sex characteristics of males. Ultrasound exam had detected no uterus and definite ovarian structure, but narrow band vaginal hypoecho and curved cavernoid structure. The child was found to have a 46,XY karyotype without an AZF deletion. DNA sequencing revealed that she has harbored a maternally derived c.323delA (p.Q108Rfs*188) variant in the nuclear receptor subfamily 5 group A member 1 (NR5A1) gene, which may result in a truncated protein. The variant was classified as pathogenic (PVS1+PM2_Supporting+PP4) based on the guidelines from the American College of Medical Genetics and Genomics. CONCLUSION: The NR5A1: c.323delA variant probably underlay the pathogenesis of 46,XY DSD in this child. The discovery of the novel variant has enriched the mutational spectrum of the NR5A1 gene and provided a basis for clinical diagnosis, treatment and prenatal diagnosis.

Clinical and genetic characteristics of a large international cohort of individuals with rare NR5A1/SF-1 variants of sex development.

BACKGROUND: Steroidogenic factor 1 (SF-1/NR5A1) is essential for human sex development. Heterozygous NR5A1/SF-1 variants manifest with a broad range of phenotypes of differences of sex development (DSD), which remain unexplained. METHODS: We conducted a retrospective analysis on the so far largest international cohort of individuals with NR5A1/SF-1 variants, identified through the I-DSD registry and a research network. FINDINGS: Among 197 individuals with NR5A1/SF-1 variants, we confirmed diverse phenotypes. Over 70% of 46, XY individuals had a severe DSD phenotype, while 90% of 46, XX individuals had female-typical sex development. Close to 100 different novel and known NR5A1/SF-1 variants were identified, without specific hot spots. Additionally, likely disease-associated variants in other genes were reported in 32 individuals out of 128 tested (25%), particularly in those with severe or opposite sex DSD phenotypes. Interestingly, 48% of these variants were found in known DSD or SF-1 interacting genes, but no frequent gene-clusters were identified. Sex registration at birth varied, with <10% undergoing reassignment. Gonadectomy was performed in 30% and genital surgery in 58%. Associated organ anomalies were observed in 27% of individuals with a DSD, mainly concerning the spleen. Intrafamilial phenotypes also varied considerably. INTERPRETATION: The observed phenotypic variability in individuals and families with NR5A1/SF-1 variants is large and remains unpredictable. It may often not be solely explained by the monogenic pathogenicity of the NR5A1/SF-1 variants but is likely influenced by additional genetic variants and as-yet-unknown factors. FUNDING: Swiss National Science Foundation (320030-197725) and Boveri Foundation Zürich, Switzerland.

Pituitary neuroendocrine tumors with PIT1/SF1 co-expression show distinct clinicopathological and molecular features.

Pituitary neuroendocrine tumors (PitNETs) are classified according to cell lineage, which requires immunohistochemistry for adenohypophyseal hormones and the transcription factors (TFs) PIT1, SF1, and TPIT. According to the current WHO 2022 classification, PitNETs with co-expression of multiple TFs are termed "plurihormonal". Previously, PIT1/SF1 co-expression was prevailingly reported in PitNETs, which otherwise correspond to the somatotroph lineage. However, little is known about such tumors and the WHO classification has not recognized their significance. We compiled an in-house case series of 100 tumors, previously diagnosed as somatotroph PitNETs. Following TF staining, histopathological features associated with PIT1/SF1 co-expression were assessed. Integration of in-house and publicly available sample data allowed for a meta-analysis of SF1-associated clinicopathological and molecular features across a total of 270 somatotroph PitNETs. The majority (74%, 52/70) of our densely granulated somatotroph PitNETs (DGST) unequivocally co-expressed PIT1 and SF1 (DGST-PIT1/SF1). None (0%, 0/30) of our sparsely granulated somatotroph PitNETs (SGST) stained positive for SF1 (SGST-PIT1). Among DGST, PIT1/SF1 co-expression was significantly associated with scarce FSH/LH expression and fewer fibrous bodies compared to DGST-PIT1. Integrated molecular analyses including publicly available samples confirmed that DGST-PIT1/SF1, DGST-PIT1 and SGST-PIT1 represent distinct tumor subtypes. Clinicopathological meta-analyses indicated that DGST-PIT1 respond more favorably towards treatment with somatostatin analogs compared to DGST-PIT1/SF1, while both these subtypes show an overall less aggressive clinical course than SGST-PIT1. In this study, we spotlight that DGST with co-expression of PIT1 and SF1 represent a common, yet underrecognized, distinct PitNET subtype. Our study questions the rationale of generally classifying such tumors as "plurihormonal", and calls for a refinement of the WHO classification. We propose the term "somatogonadotroph PitNET".

FOXL2 and NR5A1 induce human fibroblasts into steroidogenic ovarian granulosa-like cells.

Human granulosa cells in different stages are essential for maintaining normal ovarian function, and granulosa cell defect is the main cause of ovarian dysfunction. To address this problem, it is necessary to induce functional granulosa cells at different stages in vitro. In this study, we established a reprogramming method to induce early- and late-stage granulosa cells with different steroidogenic abilities. We used an AMH-fluorescence-reporter system to screen candidate factors for cellular reprogramming and generated human induced granulosa-like cells (hiGC) by overexpressing FOXL2 and NR5A1. AMH-EGFP+ hiGC resembled human cumulus cells in transcriptome profiling and secreted high levels of oestrogen and progesterone, similar to late-stage granulosa cells at antral or preovulatory stage. Moreover, we identified CD55 as a cell surface marker that can be used to isolate early-stage granulosa cells. CD55+ AMH-EGFP- hiGC secreted high levels of oestrogen but low levels of progesterone, and their transcriptome profiles were more similar to early-stage granulosa cells. More importantly, CD55+ hiGC transplantation alleviated polycystic ovary syndrome (PCOS) in a mouse model. Therefore, hiGC provides a cellular model to study the developmental program of human granulosa cells and has potential to treat PCOS.

The genetic spectrum of a Chinese series of patients with 46, XY disorders of the sex development.

PURPOSE: The etiology of 46, XY disorders of sex development (46, XY DSD) is complex, and studies have shown that different series of patients with 46, XY DSD has different genetic spectrum. In this study, we aimed to investigate the underlying genetic etiology in a Chinese series of patients with 46, XY DSD by whole exome sequencing (WES). METHODS: Seventy patients with 46, XY DSD were enrolled from the Peking Union Medical College Hospital (Beijing, China). The detailed clinical characteristics were evaluated, and peripheral blood was collected for WES to find the patients' rare variants (RVs) of genes related to 46, XY DSD. The clinical significance of the RVs was annotated according to American College of Medical Genetics and Genomics (ACMG) guidelines. RESULTS: A total of 57 RVs from nine genes were identified in 56 patients with 46, XY DSD, which include 21 novel RVs and 36 recurrent RVs. Based on the American ACMG guidelines, 43 variants were classified as pathogenic(P) or likely pathogenic (LP) variants and 14 variants were defined as variants of uncertain significance (VUS). P or LP variants were identified in 64.3% (45/70) patients of the series. Thirty-nine, 14, and 4 RVs were involved in the process of androgen synthesis and action, testicular determination and developmental process, and syndromic 46, XY DSD, respectively. The top three genes most frequently affected to cause 46, XY DSD were AR, SRD5A2, and NR5A1. Seven patients were found harboring RVs of the 46, XY DSD pathogenic genes identified in recent years, namely DHX37 in four patients, MYRF in two patients, and PPP2R3C in one patient. CONCLUSION: We identified 21 novel RVs of nine genes, which extended the genetic spectrum of 46, XY DSD pathogenic variants. Our study showed that 60% of the patients were caused by AR, SRD5A2 or NR5A1 P/LP variants. Therefore, polymerase chain reaction (PCR) amplification and Sanger sequencing of these three genes could be performed first to identify the pathogeny of the patients. For those patients whose pathogenic variants had not been found, whole-exome sequencing could be helpful in determining the etiology.

Bile acids regulate SF-1 to alter cholesterol balance in adrenocortical cells via S1PR2.

Glucocorticoid synthesis typically occurs in adrenal cortex and is influenced by cholesterol balance, since cholesterol is the sole precursor of steroids. Bile acids as the signaling molecules, have been shown to promote steroidogenesis in steroidogenic cells. However, whether bile acids directly regulate cholesterol balance remains elusive. In this study, we prepared cholestatic mouse models and cultured human adrenocortical cells (H295R) treated with taurochenodeoxycholic acid (TCDCA) to determine transcription levels of cholesterol metabolism associated genes and cholesterol concentrations in adrenocortical cells. Results showed that common bile duct ligation (CBDL) and chenodeoxycholic acid (CDCA) feeding elevated the mRNA levels of Abca1, Cyp51, Hmgcr, Srb1, and Mc2r in adrenals of mice. Meanwhile, the concentrations of total cholesterol and cholesteryl ester in adrenals of CBDL and CDCA-fed mice were dramatically lowered. The total and phosphorylation levels of HSL in adrenal glands of CBDL mice were also enhanced. Similarly, TCDCA treatment in H295R cells decreased intracellular concentrations of total cholesterol and cholesteryl ester and increased transcription levels of SRB1, MC2R, and HSL as well. Inhibition of bile acids' receptor sphingosine 1-phosphate receptor 2 (S1PR2), extracellular signal-regulated kinase (ERK) phosphorylation, and steroidogenic factor 1 (SF-1) respectively successfully abolished effect of TCDCA on H295R cells. SF-1s was found to be phosphorylated at Thr75 in TCDCA-treated H295R cells. While a mild increase of intracellular cAMP concentration was detected upon TCDCA treatment, inhibition of PKA activity with Rp-Isomer in H295R cells failed to decrease the expression of SF-1 and its target genes. Our findings suggest that conjugated bile acids affect cholesterol balance through regulation of SF-1 in adrenocortical cells so as to provide an adequate cholesterol supply for glucocorticoid synthesis, which improves and enriches our understanding of the mechanism whereby bile acids regulate cholesterol balance to affect adrenal function.

Steroidogenic Factor-1 form and function: From phospholipids to physiology.

Steroidogenic Factor-1 (SF-1, NR5A1) is a member of the nuclear receptor superfamily of ligand-regulated transcription factors, consisting of a DNA-binding domain (DBD) connected to a transcriptional regulatory ligand binding domain (LBD) via an unstructured hinge domain. SF-1 is a master regulator of development and adult function along the hypothalamic pituitary adrenal and gonadal axes, with strong pathophysiological association with endometriosis and adrenocortical carcinoma. SF-1 was shown to bind and be regulated by phospholipids, one of the most interesting aspects of SF-1 regulation is the manner in which SF-1 interacts with phospholipids: SF-1 buries the phospholipid acyl chains deep in the hydrophobic core of the SF-1 protein, while the lipid headgroups remain solvent-exposed on the exterior of the SF-1 protein surface. Here, we have reviewed several aspects of SF-1 structure, function and physiology, touching on other transcription factors that help regulate SF-1 target genes, non-canonical functions of SF-1, the DNA-binding properties of SF-1, the use of mass spectrometry to identify lipids that associate with SF-1, how protein phosphorylation regulates SF-1 and the structural biology of the phospholipid-ligand binding domain. Together this review summarizes the form and function of Steroidogenic Factor-1 in physiology and in human disease, with particular emphasis on adrenal cancer.

Correlation between LHCGR and NR5A1 genes polymorphism and male infertility risk.

INTRODUCTION: Infertility is one of the important phenomena in human reproduction. Genetic factors are the most important cause of male infertility. Here, we aimed to investigate the correlation between idiopathic male infertility and SNPs of the LHCGR (rs2293275) and NR5A1 (rs1057517779) genes in the Iranian-Azeri population. METHODS: This case-control study consisted of 100 males with infertility and 100 healthy males from the Iranian Azeri population. Genomic DNA isolation from whole blood samples and Tetra-primer amplification refractory mutation system-polymerase chain reaction (Tetra-ARMS-PCR) method was used for genotyping. The data analysis was performed by chi-square (χ2) and Fisher's exact tests. RESULTS: Genotyping analysis for LHCGR (rs2293275) polymorphism indicated that the frequency of C in the case group was significantly higher than in the control group (P < .05). Moreover, genotyping analysis for NR5A1 (rs1057517779) polymorphism indicated that the frequencies of the A allele and heterozygote GA genotype in the case group were significantly higher than those in the control group (P < .05). CONCLUSION: Our study demonstrated that the SNPs of LHCGR (rs2293275) and NR5A1 (rs1057517779) genes may play a critical role in male infertility in the Iranian Azeri population. However, further studies on other ethnic origins with larger sample sizes are essential for accessing more accurate results. Moreover, functional experiments might be needed to understand the role of these polymorphisms in the molecular pathways involved in male fertility.

Spleen function is reduced in individuals with NR5A1 variants with or without a difference of sex development: a cross-sectional study.

OBJECTIVE: NR5A1 is a key regulator of sex differentiation and has been implicated in spleen development through transcription activation of TLX1. Concerns exist about hypo- or asplenism in individuals who have a difference of sex development (DSD) due to an NR5A1 disease-causing variant. We aimed to assess spleen anatomy and function in a clinical cohort of such individuals and in their asymptomatic family member carriers. DESIGN: Cross-sectional assessment in 22 patients with a DSD or primary ovarian insufficiency and 5 asymptomatic carriers from 18 families, harboring 14 different NR5A1 variants. METHODS: Spleen anatomy was assessed by ultrasound, spleen function by peripheral blood cell count, white blood cell differentiation, percentage of nonswitched memory B cells, specific pneumococcal antibody response, % pitted red blood cells, and Howell-Jolly bodies. RESULTS: Patients and asymptomatic heterozygous individuals had significantly decreased nonswitched memory B cells compared to healthy controls, but higher than asplenic patients. Thrombocytosis and spleen hypoplasia were present in 50% of heterozygous individuals. Four out of 5 individuals homozygous for the previously described p.(Arg103Gln) variant had asplenia. CONCLUSIONS: Individuals harboring a heterozygous NR5A1 variant that may cause DSD have a considerable risk for functional hyposplenism, irrespective of their gonadal phenotype. Splenic function should be assessed in these individuals, and if affected or unknown, prophylaxis is recommended to prevent invasive encapsulated bacterial infections. The splenic phenotype associated with NR5A1 variants is more severe in homozygous individuals and is, at least for the p.(Arg103Gln) variant, associated with asplenism.

Long-read genome sequencing reveals a novel intronic retroelement insertion in NR5A1 associated with 46,XY differences of sexual development.

Despite significant advancements in rare genetic disease diagnostics, many patients with rare genetic disease remain without a molecular diagnosis. Novel tools and methods are needed to improve the detection of disease-associated variants and understand the genetic basis of many rare diseases. Long-read genome sequencing provides improved sequencing in highly repetitive, homologous, and low-complexity regions, and improved assessment of structural variation and complex genomic rearrangements compared to short-read genome sequencing. As such, it is a promising method to explore overlooked genetic variants in rare diseases with a high suspicion of a genetic basis. We therefore applied PacBio HiFi sequencing in a large multi-generational family presenting with autosomal dominant 46,XY differences of sexual development (DSD), for whom extensive molecular testing over multiple decades had failed to identify a molecular diagnosis. This revealed a rare SINE-VNTR-Alu retroelement insertion in intron 4 of NR5A1, a gene in which loss-of-function variants are an established cause of 46,XY DSD. The insertion segregated among affected family members and was associated with loss-of-expression of alleles in cis, demonstrating a functional impact on NR5A1. This case highlights the power of long-read genome sequencing to detect genomic variants that have previously been intractable to detection by standard short-read genomic testing.

A reappraisal of transcriptional regulation by NR5A1 and beta-catenin in adrenocortical carcinoma.

Background: Overexpression of the transcription factor NR5A1 and constitutive activation of canonical Wnt signalling leading to nuclear translocation of beta-catenin are hallmarks of malignancy in adrenocortical carcinoma (ACC). Based on the analysis of genomic profiles in H295R ACC cells, Mohan et al. (Cancer Res. 2023; 83: 2123-2141) recently suggested that a major determinant driving proliferation and differentiation in malignant ACC is the interaction of NR5A1 and beta-catenin on chromatin to regulate gene expression. Methods: I reanalyzed the same set of data generated by Mohan et al. and other published data of knockdown-validated NR5A1 and beta-catenin target genes. Results: Beta-catenin is mainly found in association to canonical T cell factor/lymphoid enhancer factor (TCF/LEF) motifs in genomic DNA. NR5A1 and beta-catenin regulate distinct target gene sets in ACC cells. Conclusion: Overall, my analysis suggests a model where NR5A1 overexpression and beta-catenin activation principally act independently, rather than functionally interacting, to drive ACC malignancy.

CRISPR/Cas9-mediated activation of NR5A1 steers female human embryonic stem cell-derived bipotential gonadal-like cells towards a steroidogenic cell fate.

The nuclear receptor subfamily 5 group A member 1 (NR5A1), encoding steroidogenic factor 1 (SF-1), has been identified as a critical factor in gonadal development in animal studies. A previous study of ours suggested that upregulation of NR5A1 during early gonadal differentiation in male (46,XY) human pluripotent stem cells steers the cells into a more mature gonadal cell type. However, the detailed role of NR5A1 in female gonadal differentiation has yet to be determined. In this study, by combining the processes of gonadal differentiation and conditional gene activation, we show that NR5A1 induction predominantly upregulates the female gonadal marker inhibin subunit α (INHA) and steroidogenic markers steroidogenic acute regulatory protein (STAR), cytochrome P450 family 11 subfamily A member 1 (CYP11A1), cytochrome P450 family 17 subfamily A member 1 (CYP17A1), hydroxy-delta-5-steroid dehydrogenase (HSD3B2) and hydroxysteroid 17-beta dehydrogenase 1 (HSD17B1). In contrast, NR5A1 induction did not seem to affect the bipotential gonadal markers gata binding protein 4 (GATA4) and Wilms' tumour suppressor 1 (WT1) nor the female gonadal markers r-spondin 1 (RSPO1) and wnt family member 4 (WNT4). Differentially expressed genes were highly associated with adrenal and ovarian steroidogenesis pathways. Moreover, time-series analysis revealed different dynamic changes between male and female induced samples, where continuously upregulated genes in female gonadal differentiation were mostly associated with adrenal steroidogenesis. Thus, in contrast to male gonadal differentiation, NR5A1 is necessary but not sufficient to steer human embryonic stem cell (hESC)-derived bipotential gonadal-like cells towards a more mature somatic, female cell fate. Instead, it seems to direct bipotential gonadal-like cells more towards a steroidogenic-like cell population. The information obtained in this study helps in elucidating the role of NR5A1 in gonadal differentiation of a female stem cell line.

Case Report: Longitudinal follow-up and testicular sperm extraction in a patient with a pathogenic NR5A1 (SF-1) frameshift variant: p.(Phe70Serfs*5).

Background: Steroidogenic factor 1 (SF-1), encoded by the nuclear receptor subfamily 5 group A member 1 (NR5A1) gene, is a transcriptional factor crucial for adrenal and gonadal organogenesis. Pathogenic variants of NR5A1 are responsible for a wide spectrum of phenotypes with autosomal dominant inheritance including disorders of sex development and oligospermia-azoospermia in 46,XY adults. Preservation of fertility remains challenging in these patients. Objective: The aim was to offer fertility preservation at the end of puberty in an NR5A1 mutated patient. Case report: The patient was born of non-consanguineous parents, with a disorder of sex development, a small genital bud, perineal hypospadias, and gonads in the left labioscrotal fold and the right inguinal region. Neither uterus nor vagina was detected. The karyotype was 46,XY. Anti-Müllerian hormone (AMH) and testosterone levels were low, indicating testicular dysgenesis. The child was raised as a boy. At 9 years old, he presented with precocious puberty treated by triptorelin. At puberty, follicle-stimulating hormone (FSH), luteinising hormone (LH), and testosterone levels increased, whereas AMH, inhibin B, and testicular volume were low, suggesting an impaired Sertoli cell function and a partially preserved Leydig cell function. A genetic study performed at almost 15 years old identified the new frameshift variant NM_004959.5: c.207del p.(Phe70Serfs*5) at a heterozygous state. He was thus addressed for fertility preservation. No sperm cells could be retrieved from three semen collections between the ages of 16 years 4 months and 16 years 10 months. A conventional bilateral testicular biopsy and testicular sperm extraction were performed at 17 years 10 months of age, but no sperm cells were found. Histological analysis revealed an aspect of mosaicism with seminiferous tubules that were either atrophic, with Sertoli cells only, or presenting an arrest of spermatogenesis at the spermatocyte stage. Conclusion: We report a case with a new NR5A1 variant. The fertility preservation protocol proposed at the end of puberty did not allow any sperm retrieval for future parenthood.

Induced pluripotent stem cell line generated from a patient with differences in sex development (DSD) and multiple genetic variants including a large deletion in NR5A1.

The nuclear receptor subfamily 5, Group A, Member 1 (NR5A1) gene encodes steroidogenic factor 1 (SF1), which is necessary for development of steroid hormone-producing tissues including the gonad and adrenal gland. An induced pluripotent stem cell line (iPSC) LCHi002-B was generated from a participant with differences (disorders) of sex development (DSD) and multiple genetic variants including a large deletion in NR5A1, and three single nucleotide changes in DYNC2H1, PDE4D, and ZFPM2. The line presented typical morphology, expressed stem cell markers, differentiated into three germ layers, had normal karyotype, was mycoplasma-free, and carried mutations in NR5A1, DYNC2H1, PDE4D, and ZFPM2.

Steroidogenic factor 1 (SF-1; Nr5a1) regulates the formation of the ovarian reserve.

The ovarian follicle reserve, formed pre- or perinatally, comprises all oocytes for lifetime reproduction. Depletion of this reserve results in infertility. Steroidogenic factor 1 (SF-1; Nr5a1) and liver receptor homolog 1 (LRH-1; Nr5a2) are two orphan nuclear receptors that regulate adult endocrine function, but their role in follicle formation is unknown. We developed models of conditional depletion of SF-1 or LRH-1 from prenatal ovaries. Depletion of SF-1, but not LRH-1, resulted in dramatically smaller ovaries and fewer primordial follicles. This was mediated by increased oocyte death, resulting from increased ovarian inflammation and increased Notch signaling. Major dysregulated genes were Iroquois homeobox 3 and 5 and their downstream targets involved in the establishment of the ovarian laminin matrix and oocyte-granulosa cell gap junctions. Disruptions of these pathways resulted in follicles with impaired basement membrane formation and compromised oocyte-granulosa communication networks, believed to render them more prone to atresia. This study identifies SF-1 as a key regulator of the formation of the ovarian reserve.

Genetic reanalysis of patients with a difference of sex development carrying the NR5A1/SF-1 variant p.Gly146Ala has discovered other likely disease-causing variations.

NR5A1/SF-1 (Steroidogenic factor-1) variants may cause mild to severe differences of sex development (DSD) or may be found in healthy carriers. The NR5A1/SF-1 c.437G>C/p.Gly146Ala variant is common in individuals with a DSD and has been suggested to act as a susceptibility factor for adrenal disease or cryptorchidism. Since the allele frequency is high in the general population, and the functional testing of the p.Gly146Ala variant revealed inconclusive results, the disease-causing effect of this variant has been questioned. However, a role as a disease modifier is still possible given that oligogenic inheritance has been described in patients with NR5A1/SF-1 variants. Therefore, we performed next generation sequencing (NGS) in 13 DSD individuals harboring the NR5A1/SF-1 p.Gly146Ala variant to search for other DSD-causing variants and clarify the function of this variant for the phenotype of the carriers. Panel and whole-exome sequencing was performed, and data were analyzed with a filtering algorithm for detecting variants in NR5A1- and DSD-related genes. The phenotype of the studied individuals ranged from scrotal hypospadias and ambiguous genitalia in 46,XY DSD to opposite sex in both 46,XY and 46,XX. In nine subjects we identified either a clearly pathogenic DSD gene variant (e.g. in AR) or one to four potentially deleterious variants that likely explain the observed phenotype alone (e.g. in FGFR3, CHD7). Our study shows that most individuals carrying the NR5A1/SF-1 p.Gly146Ala variant, harbor at least one other deleterious gene variant which can explain the DSD phenotype. This finding confirms that the NR5A1/SF-1 p.Gly146Ala variant may not contribute to the pathogenesis of DSD and qualifies as a benign polymorphism. Thus, individuals, in whom the NR5A1/SF-1 p.Gly146Ala gene variant has been identified as the underlying genetic cause for their DSD in the past, should be re-evaluated with a NGS method to reveal the real genetic diagnosis.